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Abstract 

We introduce a three independent functions variational formalism 
for stationary and non-stationary barotropic flows. This is less than the 
four variables which appear in the standard equations of fluid dynamics 
which are the velocity field v and the density p. It will be shown how 
in terms of our new variable the Euler and continuity equations can be 
integrated in the stationary case. 
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1 Introduction 

Variational principles for non-magnetic barotropic fluid dynamics an Eu- 
lerian variational principle are well known. Initial attempts to formulate 
Eulerian fluid dynamics in terms of a variational principle, were described 
by Herivel PP, Serrin Lin [^j. However, the variational principles devel- 
oped by the above authors were very cumbersome containing quite a few 
"Lagrange multipliers" and "potentials". The range of the total number 
of independent functions in the above formulations ranges from eleven to 
seven which exceeds by many the four functions appearing in the Eulerian 
and continuity equations of a barotropic flow. And therefore did not have 
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any practical use or applications. Seliger &: Whitham |^ have developed a 
variational formalism which can be shown to depend on only four variables 
for barotropic flow. Lynden-Bell &: Katz have described a variational 
principle in terms of two functions the load A (to be described below) and 
density p. However, their formalism contains an implicit definition for the 
velocity v such that one is required to solve a partial differential equation in 
order to obtain both v in terms of p and A as well as its variations. Much 
the same criticism holds for their general variational for non-barotropic flows 
in]. In this paper we overcome this limitation by paying the price of adding 
an additional single function. Our formalism will allow arbitrary variations 
and the definition of v will be explicit. Furthermore, we will show that for 
stationary flows using somewhat different three variational variables the Eu- 
ler and continuity equations may be reduced to a single non-linear algebraic 
equation for the density and thus can be integrated. 

We anticipate applications for this study both for stability analysis of 
known fluid dynamics configurations and for designing efficient numerical 
schemes for integrating the equations of fluid dynamics [71 |H1 El 1^ • 

The plan of this paper is as follows: We will review the basic equations of 
Eulerian fluid dynamics and give a somewhat different derivation of Seliger 
& Whitham's variational principle. Then we will describe the three function 
variational principle for non-stationary fluid dynamics. Finally we will give 
a different variational principle for stationary fluid dynamics and we will 
show how in terms the new variational variables the stationary equations 
of fluid dynamics that is the stationary Euler and continuity equations can 
be reduced to a solution of a single non-linear algebraic equation for the 
density. 

2 Variational principle of non-stationary fluid dy- 
namics 

Barotropic Eulerian fluids can be described in terms of four functions the 
velocity v and density p. Those functions need to satisfy the continuity and 
Euler equations: 





Vp{p) 




p 
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In which the pressure p{p) is assumed to be a given function of the density. 
Taking the curl of equation @ will lead to: 



1^ = V X (f; X J) (3) 

in which: 

S = V y. V (4) 

is the vorticity. Equation ^ describes the fact that the vorticity lines are 
"frozen" within the Eulerian flow^. 

A very simple variational principle for non-stationary fluid dynamics 
was described by Seliger & Whitham ^4 and is brought here mainly for 
completeness using a slightly different derivation than the one appearing in 
the original paper. This will serve as a starting point for the next section in 
which we will show how the variational principle can be simplified further. 
Consider the action: 

A = I Cifxdt 



C = C1 + C2 

A ^ P{lv'-e{p)), C2^u[^ + V-{pv)]-paf (5) 

in which e{p) is the specific internal energy. Obviously i', a are Lagrange 
multipliers which were inserted in such a way that the variational principle 
will yield the following equations: 

| + v.M = o 

Provided p is not null those are just the continuity equation and the 
conditions that /? is comoving. Let us take an arbitrary variational derivative 
of the above action with respect to v, this will result in: 

di^A = J d^xdtp5v -[v - Vv - aV/3] 

+ f dS ■ 5vpu (7) 



^The most general vortical flux and mass preserving flows that may be attributed to 
vortex lines were found in 
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Provided that the above boundary term vanishes, as in the case of astro- 
physical flows for which p = on the free flow boundary, or the case in 
which the fluid is contained in a vessel which induces a no flux boundary 
condition 6v ■ n = (re is a unit vector normal to the boundary), v must 
have the following form: 

jj=^=aV(3 + Vi^ (8) 

this is nothing but Clebsch representation of the flow field (see for example 
jl4j , |151 page 248] ) . Let us now take the variational derivative with respect 
to the density p, we obtain: 

/I diy 
(fxdt5p[-v^ ~ ^ ~ ~ ^ " ^^"^ 

+ j>dS-v5pv + j d^xv5p\\l (9) 

in which w = is the specific enthalpy. Hence provided that 5p vanishes 
on the boundary of the domain and in initial and final times the following 
equation must be satisfied: 

du 1 r, , ^ 

- = -.^-» (10) 

Finally we have to calculate the variation with respect to (3 this will lead us 
to the following results: 

5aA = f d^xdtSPi^^^ + V • (pav)] 
J at 

- jdS- vpa5f5 - J d^xpa6(3\tl (11) 

Hence choosing 5/3 in such a way that the temporal and spatial boundary 
terms vanish in the above integral will lead to the equation: 

^ + V. (12) 

Using the continuity equation ^ this will lead to the equation: 

Hence for p ^ both a and /? are comoving coordinates. Since the vorticity 
can be easily calculated from equation (jH)) to be: 

uj = V X v = Va xV/3 (14) 

Calculating ^ in which lo is given by equation (|14() and taking into account 
both equation (|13|) and equation © will yield equation Q. 
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2.1 Euler's equations 

We shall now show that a velocity field given by equation ^ , such that the 
functions a, (3, v satisfy the corresponding equations H6I1UI13() must satisfy 
Euler's equations. Let us calculate the material derivative of v: 

dv dVu da^„ dVB 

— = —— + — Vp + a— r- 15 

dt dt dt dt ^ ^ 

It can be easily shown that: 

dVu --dv -L du ^,1^9 X ^ dh' 

V— - VvkTT- = V(n^ -w) - Vvk- 



dt dt dxk 2 dxk 

dVP _ y^/^ ^/^ _ ng) 

dt dt dxk dxk 

In which is a Cartesian coordinate and a summation convention is as- 
sumed. Inserting the result from equations (|16j) into equation (|15j) yields: 

dv ^ , dv dd . -t.l _ff , 

-dt = + W + 

= -VvkVk + V{\v^ -w) = -^ (17) 
2 p 

This proves that the Euler equations can be derived from the action given 
in equation © and hence all the equations of fluid dynamics can be derived 
from the above action without restricting the variations in any way. Taking 
the curl of equation ((T7I) will lead to equation Q . 

2.2 Simplified action 

The reader of this paper might argue that the authors have introduced 
unnecessary complications to the theory of fluid dynamics by adding three 
more functions a, /?, v to the standard set v, p. In the following we will show 
that this is not so and the action given in equation © in a form suitable 
for a pedagogic presentation can indeed be simplified. It is easy to show 
that the Lagrangian density appearing in equation (jSJ can be written in the 
form: 

,dv dl3 , 1 ..^ ^,9 ;.2, 

^ = -p[^ + (^^ + ^{p)] + l^p[{v-vf -v\ 

+ ^ + V.(.p.l (18) 
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In which iT is a shorthand notation for Vi^ + aV/3 (see equation ©)• Thus 
£ has three contributions: 

^ £ + C^v ~\~ C-boundary 




The only term containing v is it can easily be seen that this term will 
lead, after we nullify the variational derivative, to equation © but will oth- 
erwise have no contribution to other variational derivatives. Notice that the 
term Cboundary Contains only complete partial derivatives and thus can not 
contribute to the equations although it can change the boundary conditions. 
Hence we see that equations © , equation (|10j) and equation ([T3|) can be de- 
rived using the Lagrangian density C in which v replaces v in the relevant 
equations. Furthermore, after integrating the four equations (|6I10I13|) we 
can insert the potentials a, /3, z/ into equation to obtain the physical ve- 
locity V. Hence, the general barotropic fluid dynamics problem is changed 
such that instead of solving the four equations (|H2j) we need to solve an 
alternative set which can be derived from the Lagrangian density C. 

2.3 The inverse problem 

In the previous subsection we have shown that given a set of functions 
a, P, V satisfying the set of equations described in the previous subsections, 
one can insert those functions into equation Q and equation ((Tl)) to obtain 
the physical velocity v and vorticity cD. In this subsection we will address 
the inverse problem that is, suppose we are given the quantities v and p how 
can one calculate the potentials a, /3, vl The treatment in this section will 
follow closely (with minor changes) the discussion given by Lynden-Bell & 
Katz '5 and is given here for completeness. 

Consider a thin tube surrounding a vortex line as described in figure Q 
the vorticity flux contained within the tube which is equal to the circulation 
around the tube is: 




(20) 



and the mass contained with the tube is: 
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Figure 1: A thin tube surrounding a vortex line 



in which dl is a length element along the tube. Since the vortex lines move 
with the flow by virtue of equation Q both the quantities A<^> and AM are 
conserved and since the tube is thin we may define the conserved load: 



in which the above integral is performed along the field line. Obviously the 
parts of the line which go out of the flow to regions in which p = has a null 
contribution to the integral. Since A is conserved is satisfies the equation: 



By construction surfaces of constant load move with the flow and contain 
vortex lines. Hence the gradient to such surfaces must be orthogonal to the 
field line: 



Now consider an arbitrary comoving point on the vortex line and donate 
it by i, and consider an additional comoving point on the vortex line and 
donate it by r. The integral: 




(22) 



VA • cD = 



(24) 




(25) 
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Figure 2: Surfaces of constant load 



is also a conserved quantity which we may denote following Lynden-Bell &: 
Katz 12] as the generalized metage. ^{i) is an arbitrary number which can 
be chosen differently for each vortex line. By construction: 



'dt 



0. 



(26) 



Also it is easy to see that by differentiating along the vortex line we obtain: 



V/x ■ uj = p 



(27) 



At this point we have two comoving coordinates of flow, namely A, /i obvi- 
ously in a three dimensional flow we also have a third coordinate. However, 
before defining the third coordinate we will find it useful to work not directly 
with A but with a function of A. Now consider the vortical flux ^{X) within 
a surface of constant load as described in figure [21 (the figure was given by 
Lynden-Bell &: Katz [S] ) . The flux is a conserved quantity and depends only 
on the load A of the surrounding surface. Now we define the quantity: 



a 



$(A) _ C7(A) 



27r 



2tt 



(28) 
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C(A) is the circulation along lines on this surface. Obviously a satisfies the 
equations: 

— = 0, • Va = (29) 

at 

Let us now define an additional comoving coordinate /3* since V/i is not 
orthogonal to the Co lines we can choose V/3* to be orthogonal to the Co lines 
and not be in the direction of the Va lines, that is we choose /?* not to 
depend only on a. Since both V/3* and Va are orthogonal to cJ, u; must 
take the form: 

u) = AVa X V/3* (30) 
However, using equation Q) we have: 

V -0 = VA - {VaxV/3*) = (31) 

Which implies that A is a function of a, (3*. Now we can define a new 
comoving function /? such that: 

(5 = 1^^ A{a,p'*)dp'*, ^ = (32) 
In terms of this function we recover the representation given in equation 

uj = VaxV/3 (33) 

Hence we have shown how a, (3 can be constructed for a known {7, p. Notice 
however, that (3 is defined in a non unique way since one can redefine (3 for 
example by performing the following transformation: (3 — > [3+ f{a) in which 
/(a) is an arbitrary function. The comoving coordinates a, (3 serve as labels 
of the vortex lines. Moreover the vortical flux can be calculated as: 

$ = j u-dS = J dadf3 (34) 

Finally we can use equation © to derive the function v for any point s 
within the flow: ^ 

iy{s) = j^{v- aV(3) ■ df+ (35) 

in which i is any arbitrary point within the flow, the result will not depend 
on the trajectory taken in the case that u is single valued. If u is not single 
valued on should introduce a cut, which the integration trajectory should 
not cross. 
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2.4 Stationary fluid dynamics 

Stationary flows are a unique phenomena of Eulerian fluid dynamics which 
has no counter part in Lagrangian fluid dynamics. The stationary flow is 
defined by the fact that the physical fields v, p do not depend on the temporal 
coordinate. This however does not imply that the stationary potentials 
a, (3, V are all functions of spatial coordinates alone. Moreover, it can be 
shown that choosing the potentials in such a way will lead to erroneous 
results in the sense that the stationary equations of motion can not be 
derived from the Lagrangian density £ given in equation H19() . However, 
this problem can be amended easily as follows. Let us choose a, v to depend 
on the spatial coordinates alone. Let us choose /? such that: 

P = f3-t (36) 

in which /3 is a function of the spatial coordinates. The Lagrangian density 
C given in equation (|19)) will take the form: 

£ = p(a-e(p)-i(W + aV/3)2) (37) 

Varying the Lagrangian L = J LSx with respect to a, /3, p leads to the 
following equations: 

V • {p^}) = 
pv ■ Va = 
p(^-V/3-l) = 

1 -2 

a = -v +h (38) 

a is thus the Bernoulli constant (this was also noticed in j7j). Calcula- 
tions similar to the ones done in previous subsections will show that those 
equations lead to the stationary Euler equations: 

p{ij-V)i; = -Vp{p) (39) 

3 A simpler variational principle of non- stationary 
fluid dynamics 

Lynden-Bell &: Katz |S] have shown that an Eulerian variational principle 
for non-stationary fluid dynamics can be given in terms of two functions the 
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density p and the non-magnetic load A defined in equation (|22p . However, 
their velocity was given an implicit definition in terms of a partial differen- 
tial equation and its variations was constrained to satisfy this equation. In 
this section we will propose a three function variational principle in which 
the variations of the functions are not constrained in any way, part of our 
derivation will overlap the formalism of Lynden-Bell & Katz. The three vari- 
ables will include the density /o, the non-magnetic load A and an additional 
function to be defined in the next subsection. This variational principle is 
simpler than the Seliger & Whitham variational principle jlj which is given 
in terms of four functions and is more convenient than the Lynden-Bell &: 
Katz ^ variational principle since the variations are not constrained. 

3.1 Velocity representation 

Consider equation ()24() . since to is orthogonal to VA we can write: 

w = i? X VA (40) 

in which K is some arbitrary vector field. However, since V - a; = it follows 
that K = V0 for some scalar function theta. Hence we can write: 

Vxt; = a; = WxVA (41) 

This will lead to: 

v = dVX + Vu (42) 

For the time being v is an arbitrary scalar function, the choice of nota- 
tion will be justified later. Consider now equation (|23)1 . inserting into this 
equation v given in equation (|42|) will result in: 

^ = ^ + t;.VA = ^ + [eV\ + Vu) . VA = 0. (43) 
at at at 

This can be solved for 9, the solution obtained is: 



V ivA|2 ; 

Inserting the above expression for 9 into equation (|42() will yield: 

dx dx 



(44) 



A + Vi/ - A(A • Vzv) = — ^A + W (45) 
|VA| |VA| 
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in which A = is a unit vector perpendicular to the load surfaces and 

V*z^ = Vj^ — A(A • Vi^) is the component of Vz^ parallel to the load surfaces. 
Notice that the vector v — Vu is orthogonal to the load surfaces and that: 

\v-Vu\ = {v-Viy)-X = d\VX\^9= (^"Y'^)-^ (46) 

|VA| 

Further more by construction the velocity field v given by equation (|^^ 
ensures that the load surfaces are comoving. Let us calculate the circulation 
along A surfaces: 

C(A) = (p v-df= (t) V*u ■dr= (l Vu-df= [u]x (47) 
Jx Jx Jx 

[iy]x is the discontinuity of v across a cut which is introduced on the A 
surface. Hence in order that circulation C(A) on the load surfaces (and 
hence everywhere) will not vanish v must be multiple-valued. Following 
Lamb ^15- page 180, article 132, equation 1] we write v in the form: 

u = C7(A)P, [p]a = 1 (48) 

in terms of P the velocity is given as: 

dx 

v = — I^A + C7(A)V*z? (49) 
|VA| 

And the explicit dependence of the velocity field v on the circulation along 
the load surfaces C(A) is evident. 

3.2 The variational principle 

Consider the action: 

A = I Cd^xdt 



C = C1 + C2 

A ^ p{\v^-e{p)), C2^u[^ + V-{pv)] (50) 

In which v is defined by equation (|45|) . u is not a simple Lagrange multiplier 
since v is dependent on u through equation H45|) . Taking the variational 
derivative of C with respect to 1^ will yield: 

6^C = 6iy[^+V ■ {pv)] + pv ■ 6^v + uV • {p6^v) (51) 
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This can be rewritten as: 

5^C = 5v[^+V ■ {pv)] + p{v - Vu) -b^v + V- {pu6^v) (52) 
Now by virtue of equation : 

6^v = V*6u (53) 

which is parallel to the load surfaces, while from equation (|42l) we see that 
v — Vu is orthogonal to the load surfaces. Hence, the scalar product of those 
vector must be null and we can write: 

6^C = 6iy[^ + V • {pv)] + V • ipi^V*du) (54) 

Thus the action variation can be written as: 

dp 



S^A = J d^xdt5v[-^ + V • {pv)] 



+ jdS-pvV*5v (55) 

This will yield the continuity equation using the standard variational pro- 
cedure. Notice that the surface should include also the "cut" since the 
function is in general multi valued. Let us now take the variational derivative 
with respect to the density p, we obtain: 

5pA = j d^xdt5p\^v^ - w - ^ - V ■ Vu] 

+ j>dS-v5pv + j d^xu5p\\l (56) 

Hence provided that 5p vanishes on the boundary of the domain and in 
initial and final times the following equation must be satisfied: 

— = -V — w 57 
dt 2 ^ ' 

This is the same equation as equation (|1U() and justifies the use of the sym- 
bol V in equation H42|) . Finally we have to calculate the variation of the 
Lagrangian density with respect to A this will lead us to the following re- 
sults: 

5\C = pv ■ 6xv + z/V • {p5\v) = p{v — Vz/) • 5\v + V • {pudxv) 

= pe\V\\{\-5xv)+V -{pvSxv) (58) 
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in equation was used. Let us calculate 6\v, after some straightforward 
manipulations one arrives at the result: 



6xv 



A 



|VA| 



d{6X) 



(59) 



Inserting equation into equation (|55)) and integrating by parts will yield: 

-d{pe) 



5xC = SX 



dt 



+ V • {pOv) 



V • [p{5xv,. - 9v5X)] - ^^^^ (60) 



Hence the total variation of the action will become: 



SxA 



(TxdtSX 



d{pe) 

dt 



+ V • {pOv) 



+ jdS- {5xvu - ev5X)p - / d'xpOSXWl 



(61) 



Hence choosing 6X in such a way that the temporal and spatial boundary 
terms vanish in the above integral will lead to the equation: 



d{pe) 

dt 



+ V • {pOv) = 



Using the continuity equation (P) will lead to the equation: 



dt 







(62) 



(63) 



Hence for p 7^ both A and 9 are comoving. Comparing equation (jH)) to 
equation we see that a is analogue to 6 and /3 is analogue to A and 
all those variables are comoving. Furthermore, the v function in equation 
1)42(1 satisfies the same equation as the v appearing in equation Q which is 
equation (|57|) . It follows immediately without the need for any additional 
calculations that v given in equation (|^^ satisfies Euler's equations Q, the 
proof for this is given in subsection 12.11 in which one should replace a with 
9 and (5 with A. Thus all the equations of fiuid dynamics can be derived 
from the action ()5U() without restricting the variations in any way. The 
reader should notice an important difference between the current and pre- 
vious formalism. In the current formalism is a dependent variable defined 
by equation (|44|) . while in the previous formalism the analogue quantity a 
was an independent variational variable. Thus equation (|6c{() should be con- 
sidered as some what complicated second-order partial differential equation 
(in the temporal coordinate t) for A which should be solved simultaneously 
with equation (j^7|) and equation (P). 
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3.3 Simplified action 

The Lagrangian density £ given in equation (|5U|) can be written explicitly 
in terms of the three variational variables p,X,iy as follows: 



£, 



C 

L 

boundary 




-iVvf 



dv 
'dt 



e{p) 



(64) 



Notice that the term Choundary contains only complete partial derivatives and 
thus can not contribute to the equations although it can change the bound- 
ary conditions. Hence we see that equation (P), equation ((S7j) and equation 
can be derived using the Lagrangian density C in which v is given in 
terms of equation (jlSJ) in the relevant equations. Furthermore, after inte- 
grating those three equations we can insert the potentials A, v into equation 
H45() to obtain the physical velocity v. Hence, the general barotropic fluid 
dynamics problem is altered such that instead of solving the four equations 
H1I2|) we need to solve an alternative set of three equations which can be de- 
rived from the Lagrangian density C. Notice that the specific choice of the 
labelling of the A surfaces is not important in the above Lagrangian density 
one can replace: A— > A(A), without changing the Lagrangian functional 
form. This means that only the shape of the A surface is important not 
their labelling. In group theoretic language this implies that the Lagrangian 
is invariant under an infinite symmetry group and hence should posses an 
infinite number of constants of motion. In terms of the Lamb type function 
v' defined in equation (|48|) . the Lagrangian density given in equation (|5fl|) 
can be rewritten in the form: 



£ = 
£ = 

boundary = 



~l~ ^boundary 

1 / § + C(A)VA- VP 



dt 



|VA| J 
V • {C{X)i?pv) 



^{C{X)Vuf-CiX)^-e{p) 



(65) 



Which emphasize the dependence of the Lagrangian on the the circulations 
along the load surfaces C(A) which are given as initial conditions. 
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3.4 Stationary fluid dynamics 



For stationary flows we assume that both the density p and the load A are 
time independent. Hence the velocity field given in equation ()45|) can be 
written as: 

^ = Vz^ - A(A • Vu) = V*u = C{X)V*v (66) 

thus the stationary flow is parallel to the load surfaces. From the above 
equation we see that in the stationary case v can be written in the form: 



1^ = 1^0- /(A,t) 



(67) 



in which /(A,t) is an arbitrary function and z/q is independent of the tem- 
poral coordinate. Hence we can rewrite the velocity v as: 

iJ= VVo = C7(A)V*z^o (68) 
Inserting equation 1)6 7|) and equation H68|) into equation (jSTf) will yield: 

dfiX,t) 1 



dt 



^v'^ + w = B{X) 



(69) 



in which -B(A) is the Bernoulli constant. Integrating we obtain: 

fiX,t) = B{X)t + g{X) (70) 

the arbitrary g{X) function can be absorbed into vq and thus we rewrite 
equation 1)6 7|) in the form: 



1/ = i/Q - B{X)t 



(71) 



Further more we can rewrite the conserved quantity 9 given in equation (|44|) 
as: 



VA- Vi/ 

. |VA|2 , 



A • Vfo 

. |VA| , 



+ t 



dB{X) 
dX 



(72) 



The Lagrangian density C given in equation (|50j) can be written in the sta- 
tionary case taking into account equation (|^S)) and equation (f7T|) as follows: 



C 



i(A.w)2-i(W)2 + S(A) -£(/,) 
B{X)-^{V*uof-e{p) 
B{X)-kc{X)V*i?of-eip) 



(73) 
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Taking the variational derivative of the Lagrangian density L with respect to 
the mass density p will yield the Bernoulli equation H69|) . The variation of the 
Lagrangian L = J d^xC with respect to will yield the mass conservation 
equation: 

V* • {pv) = (74) 

this form is equivalent to the standard stationary continuity equation 

V • (pv) = since there is no mass flux orthogonal to the load surfaces. 

Finally taking the variation of L with respect to A will yield: 



dX \ IVAI 







(75) 



Which can be also obtained by inserting equation H72|) into equation H63() . 
Hence we obtained three equations (|69I74I75|) for the three spatial functions 
p,X and 1^0 . Admittedly those equations do not have a particularly simple 
form, we will obtain a somewhat better set of equations in the next section. 



4 Simplified variational principle for stationary fiuid 
dynamics 

In the previous sections we have shown that fluid dynamics can be described 
in terms of four first order differential equations and in term of an action 
principle from which those equations can be derived. An alternative deriva- 
tion in terms of three differential equations one of which (equation (|63() 1 is 
second order has been introduced as well. Those formalisms were shown to 
apply to both stationary and non-stationary fluid dynamics. In the following 
a different three functions formalism for stationary fluid dynamics is intro- 
duced. In the suggested representation, the Euler and continuity equations 
can be integrated leaving only an algebraic equation to solve. 

Consider equation (jl^j) . for a stationary flow it takes the form: 



Hence v can take the form: 



v-Va = (76) 



^ VaxK 

V = (77) 

P 



However, since the velocity field must satisfy the stationary mass conserva- 
tion equation equation (^: 

V • (pv) = (78) 
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We see that K must have the form K = VN, where N is an arbitrary 
function. Thus, v takes the form: 

v = (79) 

P 

Let us now calculate v x Cu in which uj is given by equation (fTl|) . hence: 

iTxu; = ( ) X (Va X V/3) 

P 

= -Va(Va X ViV) • V/3 (80) 

Now since the flow is stationary can be at most a function of the three 
comoving coordinates a,/3,^ defined in subsections 12.21 and 12.41 hence: 

dN ^ dN ^- dN ^ 

Inserting equation (|HT|) into equation (|5n)) will yield: 

I ^ ON ^ 

iJ xuj = -Va—(VaxVfi) -VP (82) 
P op 



Rearranging terms and using vorticity formula (|14j) we can simplify the 
above equation and obtain: 

^ ^ 1 ^ SAT ^ ^ 

vxio = — Va— (V/x-cj) (83) 
P dn 

However, using equation ()27() this will simplify to the form: 

- dN 

vXio = -Va— (84) 
oil 

Now let us consider equation (jS)), for stationary flows this will take the form: 

V X (iJ X w) = (85) 
Inserting equation (IM)) into equation (|85|) will lead to the equation: 

^ dN 

V(^) X Va = (86) 
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However, since N is at most a function of a, /?, /u. It follows that ^ is some 
function of a: 

dN 

= -^(°' '«^> 

This can be easily integrated to yield: 

N = -nF{a) + G{a,/3) (88) 

Inserting this back into equation (|79|) will yield: 

Va X (-F(a)V/i + Ifv^) 

i/= ^ (89) 

P 

Let us now replace the set of variables a,f3 with a new set a',f3' such that: 



a 



F{a)da, p' = (90) 
Fia) 



This will not have any effect on the vorticity presentation given in equation 
(dU since: 

c5 = Va X V/3 = Va X V/3 = Va X Vp' (91) 

However, the velocity will have a simpler presentation and will take the 
form: 

^ ^ Va' xV{-f, + G'ia'J')) 
P 

in which G' = ^. At this point one should remember that was defined 
in equation H25() up to an arbitrary constant which can very between vortex 
lines. Since the lines are labelled by their a',/3' values it follows that we 
can add an arbitrary function of a',/5' to /i without effecting its properties. 
Hence we can define a new /i' such that: 

= fi-G'{a',P') (93) 

Inserting equation (jHS)) into equation (|^^ will lead to a simplified equation 
for v: 

^ Va' X Va' , 

V = — 94 

P 

In the following the primes on a, (3, fi will be ignored. It is obvious that v 
satisfies the following set of equations: 

iT • V^ = 0, v-Va = 0, v-Vf3 = l (95) 
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to derive the right hand equation we have used both equation (|27|) and 
equation dJ. Hence ;U,a are both comoving and stationary. As for /3 it 
satisfies equation (|38() . 

By vector multiplying v and uj and using equations H94I14() we obtain: 

V X uj = Va (96) 

this means that both v and uj lie on a surfaces and provide a vector basis 
for this two dimensional surface. 



4.1 The action principle 

In the previous subsection we have shown that if the velocity field v is given 
by equation (|94j) than equation is satisfied automatically for stationary 
flows. To complete the set of equations we will show how the Euler equations 
((21) can be derived from the Lagrangian: 

L = j Cd^x 

C ^ p{lv'-e{p)) (97) 

In which v is given by equation ()94|1 and the density p is given by equation 
(|2Z1): 

p = V/x • d; = V^i • (Va X V/3) = (98) 

d{x,y,z) 

In this case the Lagrangian density of equation ()97() will take the form: 

£ = ,(i{5^i£l2:)^ (98) 

and can be seen explicitly to depend on only three functions. The variational 
derivative of L given in equation 1)9 7|) is: 

5L = j 5Cd-^x 

SC = 5p{v^ - w{p)) +pv-5v (100) 

Let us make arbitrary small variations 6ai = {6a, 5(3, 6fi) of the functions 
= {a,(3,fi). Let us define the vector: 

^ dr 

^ = -T^5a^ (101) 
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This will lead to the equation: 

5ai = -Vai ■ i (102) 

Making a variation of p given in equation (jUHJ with respect to will yield: 

6p = -V- (pi) (103) 

(for a proof see for example ^21) • Calculating 5v by varying equation ()94|) 
will give: 

6v = -—V + -V X (pix v) (104) 
P P 

Inserting equations (|lfl3llfl4)) into equation (|10fl|) will yield: 

5C = V -V X {p^ X v) - 6p{^i)^ + w) 

= v-Vx{p(xv)+V-{pi){^v'^ + w) (105) 
Using the well known vector identity: 

A - V X {C X A) = V ■ {{C X A) X A) + {C X A) -V X A (106) 
and the theorem of Gauss we can write now equation H10U|) in the form: 

5L = (j) dS ■ [{i X v) X V + {^if + w)i]p 

+ J d^xi- [vxuj- V(^i7^ + w)]p (107) 

Suppose now that 5L = for a such that the boundary term in the above 
equation is null but that ^ is otherwise arbitrary, then it entails the equation: 

pvxCo- pV{-v^ + w) = (108) 
Using the vector identity : 

^V{v'^) = {v-V)v + vx{V xv) (109) 
and rearranging terms we recover the stationary Euler equations: 

p{v-V)v = -Vp (110) 
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4.2 The integration of the stationary Euler equations 

Let us now combine equation (|1U8|) with equation (|96|). this will yield: 

^a = V{^v^ + w) (111) 

which can easily be integrated to yield: 

a = ^v'^ + w{p) (112) 

in which an arbitrary constant is absorbed into a, this is consistent with 
equation (|38|) and implies that a is the Bernoulli constant. Let us now 
insert v given by equation into the above equation: 

1 / V/i X VaX ^ , , , ^ 

»=2[ p ) +^('') 

This can be rewritten as: 

w{p)p'^ - p^a + ^{VpxVaf = (114) 

By defining: 

U ^ p' 

W{U) = w{p) (115) 

We obtain: 

W{U)U-Ila + ^{VpxVaf =0 (116) 

If one chooses an a and p functions such that the above nonlinear equation 
has a positive (or zero) solution for 11 at every point in space than one 
obtains a solution of both the Euler and continuity equations in which v is 
given by equation (|94|) and p = +V^. One should notice, however, that 
a and p are not arbitrary functions, a is a function of the load defined in 
equation (|22|) and thus according to equation (|24j) must satisfy the equation: 

Va-w = (117) 

that is both the velocity and vorticity fields must lie on the alpha surfaces. 
While for the velocity v this is assured by equation (|94() this is not assured 
for the vorticity. Furthermore, p is not an arbitrary function it is a metage 
function an thus must satisfy equation H27() . Moreover, obtaining a solution 
does not imply that the obtained solution is stable. Having said that we 
notice that the technique requires only the solution of the algebraic equation 
(|116|1 and does involve solving any differential equations. 
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5 Conclusion 



In this paper we have reviewed Eulerian variational principles for non- 
stationary barotropic fluid dynamics and introduced a simpler three inde- 
pendent functions variational formalisms for stationary and non-stationary 
barotropic flows. This is less than the four variables which appear in the 
standard equations of fluid dynamics which are the velocity field v and the 
density p. We have shown how in terms of our new variables the stationary 
Euler and continuity equations can be integrated, such that the station- 
ary Eulerian fluid dynamics problem is reduced to a non-linear algebraic 
equation for the density. 

The problem of stability analysis and the description of numerical schemes 
using the described variational principles exceed the scope of this paper. We 
suspect that for achieving this we will need to add additional constants of 
motion constraints to the action as was done by see also ^Sl) hope- 

fully this will be discussed in a future paper. 
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